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Abstract.—Life history studies were conducted on four commonly occurring species of 
mayflies from an Ozark stream in Missouri. A multivariate statistical technique, using a 
suite of In-transformed body measurements of Stenonema mediopunctatum (Mc- 
Dunnough) (Ephemeroptera: Heptageniidae) nymphs, discerned at least nine instars during 
the sampling period. This technique provided greater resolution in discriminating among 
instars than using raw data or head capsule widths alone. Further, S. mediopunctatum was 
determined to be univoltine in these streams, as were Caenis hilaris (Say) and C. lati- 
pennis Banks (Ephemeroptera: Caenidae). Tricorythodes curvatus Allen (Ephemeroptera: 
Leptohyphidae) apparently was multivoltine. 
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The shift in ecological research from de- 
scriptive natural history to more quantita- 
tive approaches requires that life history 
and life cycle information constitute a larg- 
er part of studies dealing with community- 
level processes (Rosenberg 1979); however, 
life histories of most species of aquatic in- 
sects are not known in detail (Richardson 
and Tarter 1976). This type of information 
is important for a better knowledge not only 
of the insects themselves, but of their dis- 
tribution and ecology. Mayflies are an im- 
portant component of aquatic ecosystems 
and the life histories must be determined in 
order to have a complete understanding of 
their functions in aquatic ecosystems (Kon- 
dratieff and Voshell 1980). Further, knowl- 
edge of the number of instars for each may- 
fly species or population, and how this in- 
formation correlates with environmental 
factors, will lead to a better understanding 
of their biology (Fink 1982, 1984), espe- 
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cially when describing life histories, life cy- 
cles, or energy budgets (Benton and Prit- 
chard 1988). 

The nymphal lifespan of temperate may- 
flies varies from 7—13 days (Gray 1981) to 
approximately 2.5 yr, depending on the spe- 
cies and geographic location of the popu- 
lation (see Brittain 1982, Sweeney et al. 
1995). Instar number, which is variable in 
mayflies (Brittain 1982, Fink 1982, 1984), 
ranges from 10 to 50, with most species in 
the range of 15-25 (see Brittain 1982, But- 
ler 1984). Further, the number of nymphal 
instars is not constant for a particular spe- 
cies, but may vary with sex and environ- 
mental conditions, such as temperature and 
nutrition, during development (Brittain 
1990). Although voltinism has been deter- 
mined for many species, the number of in- 
stars has been published for only a few. For 
example, the number of instars for Steno- 
nema canadense Walker is between 40 and 
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45 (Ide 1935), Ephemera simulans Walker 
approximately 30 (Ide 1935), Ameletus cel- 
er McDunnough 21 (Benton and Pritchard 
1988), and S. modestum (Banks) 14-15 in 
the overwintering brood (Kondratieff and 
Voshell 1980). 

Although mayflies characteristically are 
heterogeneous in growth, morphological 
development, and instar number (Fink 
1982), studies of growth and development 
often rely on body measurements. An early 
method, introduced by Ide (1935), was 
based on the change in segmentation of the 
caudal filaments. The number of instars that 
a fully grown nymph has passed through 
may be estimated by counting the number 
of segments of the caudal filaments. How- 
ever, it is difficult to obtain specimens from 
the field with intact caudal filaments, par- 
ticularly in a quantified manner. A com- 
monly used metric for mayflies and many 
other insect taxa is that of head capsule 
width (e.g., Kondratieff and Voshell 1980), 
which after direct measurement, are cate- 
gorized into size classes. The most com- 
monly used methods for mayfly instar de- 
termination are rearing, with direct obser- 
vation of instars, and the simple frequency 
method, in which instars are indirectly de- 
termined by plotting the number of individ- 
uals per size class (Fink 1984). Another 
popular method is the Janetschek method 
which requires the calculation of a moving 
average of each simple frequency size class. 
However, head capsule width, the simple 
frequency and Janetschek methods, as well 
as others including the Cassie method and 
Dyar’s Law, have been suggested as unre- 
liable (Fink 1984). In a study that examined 
the utility of six mensural characters of 
Ephemera danica Muller (Ephemeridae) for 
determining growth and voltinism, size 
classes were established based on head cap- 
sule width (Aguayo-Corraliza et al. 1991). 
However, the data were not used to approx- 
imate instar number. In that study, only 
wingpad length was shown to exhibit allo- 
metric growth; whereas, body, head, and 
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leg lengths and widths exhibited isometric 
growth. 

Among mid-temperate mayflies, univol- 
tine life cycles are the most common (Brit- 
tain 1990), whereas tropical and subtropical 
populations can deviate from this pattern 
with shorter, aseasonal life cycles (e.g., 
Nolte et al. 1996). The duration of a life 
cycle depends in part on factors that influ- 
ence growth and development in all life 
stages. Also, genetic constraints that limit 
rates of these processes may exist (Butler 
1984). Although thermal regime can have a 
major influence on voltinism of aquatic in- 
sects, life history patterns are the result of 
a complex bioenergetic interaction between 
temperature and food abundance and qual- 
ity (Anderson and Cummins 1979, Swee- 
ney 1984). 

The family Heptageniidae is widely dis- 
tributed throughout the Holarctic region and 
comprises a diverse component of benthic 
communities in Oriental (Dudgeon 1996) 
and Palearctic streams. More specifically, 
heptageniid nymphs of the genus Sterone- 
ma frequently are the most abundant ben- 
thic insects in streams and rivers of eastern 
North America (Bednarik and McCafferty 
1979). The elucidation of the complete life 
cycle and life histories of these mayflies 
and other aquatic insects in the stream com- 
munity is fundamental for a complete un- 
derstanding of community dynamics (Rich- 
ardson and Tarter 1976). This research was 
conducted to determine if using a suite of 
nymphal morphological characters with 
multivariate statistical analyses would pro- 
vide greater resolution than other common- 
ly used techniques in estimating the number 
of instars for a given mayfly species. Fur- 
ther, to examine voltinism, head capsule 
widths were taken and size-frequencies 
were calculated for four lotic mayfly spe- 
cies. 


MATERIALS AND METHODS 


The site selected for life history studies 
was the Meramec River at the University of 
Missouri Hugo Wurdack Research Farm in 
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Crawford County, Missouri. At this locali- 
ty, the Meramec River is a 5th order stream 
located in the Meramec River Basin, which 
is in the Ozark Plateau physiographic re- 
gion. The study area is within a cleared/ 
grazed land use area with a narrow, natural 
riparian zone. Riffle habitats were quanti- 
tatively sampled with a Surber sampler, 
which has a base area of 30.5 X 30.5 cm 
(1 ft?) and mesh openings of 1 mm. The 
sediment within the frame was disturbed to 
a depth of approximately 6 cm with a hand 
rake for approximately | min. Benthic in- 
vertebrates collected into the net of the 
sampler were transferred into bottles con- 
taining 80% ethyl alcohol and transported 
to the laboratory. The sampling regime for 
1992-93 consisted of weekly collections of 
40 samples from mid-May through late 
September, and monthly collections for the 
remainder of the sampling year, for a total 
of 28 sampling dates. 

Instar determination of Stenonema me- 
diopunctatum (McDunnough).—Of the 
mayfly nymphs collected, S. mediopuncta- 
tum occurred consistently and in large num- 
bers, making it a candidate taxon for anal- 
ysis of body measurements and instar de- 
termination. Nymphs of S. mediopunctatum 
from representative subsamples of 13 of the 
28 collecting dates, from May through No- 
vember 1992, were measured (mm) using 
an ocular micrometer. Measurements con- 
sisted of head length and width, body 
width, and lengths of the profemur, protibia, 
mesofemur, mesotibia, metafemur, and me- 
tatibia. Head capsule was measured at the 
greatest width, which included the com- 
pound eyes. Body width also was measured 
at its greatest width, which was the prono- 
tum. A total of 322 specimens was exam- 
ined. 

Cluster analyses were performed for head 
capsule width separately, and for the entire 
data set, which included all nine mensural 
characters. Hierarchical methods of classi- 
fication, or cluster analyses, operate on a 
matrix of similarities among a set of units. 
This technique places variables that are 
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highly correlated and similar to each other 
into groups and excludes from clusters 
those variables that are unlike. With the av- 
erage linkage method (UPGMA), the dis- 
tance between two clusters is the average 
distance between pairs of observations, one 
in each cluster. This analytical technique is 
a commonly used tool in morphometric re- 
search (e.g., Strauss 1992, Wool and Man- 
heim 1992). 

Because shape change accounts for sig- 
nificant variation among instars of insects 
undergoing allometric development, the 
data were In-transformed and all body mea- 
surements again were clustered. Only a sin- 
gle dimension can be examined with linear 
measurements; therefore, the In-trans- 
formed values were used to provide a math- 
ematical approximation of shape by calcu- 
lating two-dimensional variables (In A + In 
B = In (A * B)I. 

After cluster analysis had grouped the in- 
dividuals, first-level (most similar) group- 
ings were assumed to represent instars and 
were numbered. These instar numbers were 
then assigned to all individuals in the orig- 
inal database. Discriminant function analy- 
sis (DFA) simultaneously maximizes inter- 
group differences and minimizes intragroup 
variation among individuals by altering the 
linear combination of variables on each of 
a number of orthogonal axes. DFA then was 
used to determine the degree of overlap 
among the clustered groups (instars), giving 
an indication of the range of variation with- 
in and among instar groupings. The subse- 
quent classification phase of DFA then as- 
signs each specimen to an instar based on 
the linear combination of variables from 
each discriminant function axis. Percent of 
correct assignments was used as a separate 
measure of morphometric distinction 
among instars. More specifically, percent 
correct classification was used to determine 
the accuracy of instar assignments. All sta- 
tistical analyses were conducted using 
SPSS, version 4.0 (SPSS Inc., Chicago, IL). 

Voltinism determination of lotic may- 
flies.—The four species chosen for analysis 
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of voltinism were those with nymphs that 
were collected throughout the year: S. me- 
diopunctatum (Heptageniidae), Tricorytho- 
des curvatus Allen (Leptohyphidae), Caenis 
hilaris (Say), and C. latipennis Banks 
(Caenidae). A black light was used to col- 
lect adult mayflies, and was positioned near 
the stream 12 h prior to each collection pe- 
riod. Using an ocular micrometer, nymphal 
head capsule measurements were tabulated 
for approximately 600-700 specimens of 
each of the four species from 17 of the 28 
collecting dates (May 1992 through April 
1993). The number of specimens collected 
and measured per month for $. mediopunc- 
tatum ranged from 15 to 54, whereas T. 
curvatus ranged from 4 to 50, C. hilaris 
from 3 to 50, and C. latipennis from 6 to 
50. Size-frequency histograms then were 
constructed to determine voltinism. 


RESULTS 


Instar determination of Stenonema me- 
diopunctatum.—Using cluster analysis, raw 
data of head capsule width measurements 
clustered into five distinct groups at a small 
dissimilarity distance (Fig. 1A). The DFA 
classification rate was 89.6% using this ap- 
proach. When raw data of all mensural 
characters were clustered, the data again ex- 
hibited five distinct groups (Fig. 1B), how- 
ever the DFA classification rate improved 
to 95.6%. Finally, In-transformed data for 
all mensural characters clustered into nine 
groups (Fig. 1C), and still maintained a 
high level of DFA classification (95.3%), 
indicating very distinct groupings. 

Voltinism determination of lotic may- 
flies —Stenonema mediopunctatum showed 
mainly incosistent size-frequency distribu- 
tions throughout the year, although some 
patterns can be observed (Fig. 2A). The 
largest head capsule width measurements 
occurred in May and April for 1992 and 
1993, respectively, although the frequencies 
of these measurements were low. Numbers 
remained high throughout the growing sea- 
son for this species, and adults emerged 
during June, July, and August. 
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Tricorythodes curvatus exhibited very 
small changes in size-frequency distribu- 
tions throughout the year (Fig. 2B). From 
May through September, head capsule 
widths were as high as 1 mm; however, 
from October through April, the population 
was represented primarily by higher fre- 
quencies of smaller individuals. Adults of 
this species emerged during June, July, Au- 
gust, and September, and a second cohort 
can be observed in September. 

Caenis hilaris was the least abundant of 
the four species examined, and no speci- 
mens were collected during December, Jan- 
uary, February, or March (Fig. 2C). Peak 
abundance and the largest specimens oc- 
curred during early summer months (i.e., 
late April, May, June). Caenis latipennis 
exhibited a pattern similar to that of C. hi- 
laris, with the largest sized individuals oc- 
curring in April and May (Fig. 2D). Adult 
emergence of both species of Caenis oc- 
curred during June through September. 


DISCUSSION 


Analysis of mensural characters of S. me- 
diopunctatum nymphs revealed that at least 
nine instars were present in the samples 
from the Meramec River during the study 
period. Because a single outlier existed in 
each analysis, representing an instar mark- 
edly dissimilar from the others, it is likely 
that the morphometric distribution does not 
represent a contiguous succession of in- 
stars. Using a suite of In-transformed body 
measurements provided increased resolu- 
tion over using non-transformed data and 
increased resolution and precision over us- 
ing head capsule widths alone. Cluster anal- 
ysis for head capsule widths resulted in five 
instars and the lowest classification rate 
(89.6%) of the three methods evaluated. 
Thus, a relatively high degree of overlap 
(11.4%) existed among groups (instars) 
when based on only a single variable. When 
a suite of nine mensural variables was used, 
the number of instars remained unchanged, 
although classification accuracy increased 
dramatically to 95.6%. Thus, the number of 
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Fig. 1. Dendrogram of cluster analysis of Stenonema mediopunctatum nymphs based on measurements of 
(A) head capsule widths, (B) raw data of nine body characters, and (C) In-transformed data of nine body 
characters. Numbers on the left of each dendrogram indicate the number of individuals within a given cluster. 
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Fig. 2. Percent of individuals in each head capsule size class within each sample of four species of Ephem- 
eroptera collected from the Meramec River from March 1992 through April 1993. A) Stenonema mediopunc- 
tatum, B) Tricorythodes curvatus, C) Caenis hilaris, D) Caenis latipennis. Arrows indicate adult presence in 
black light samples, numbers in parentheses indicate the number of immatures measured. 
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misclassifications was approximately 
halved while maintaning the same level of 
grouping resolution. The use of In-trans- 
formed variables provided dramatically in- 
creased discrimination among groups, and 
cluster analysis identified nine groups. All 
characters were highly correlated with DFA 
axis | and accounted for >95% of the var- 
iation for both raw data and In-transformed 
data. Thus, all characters measured were 
important in discriminating among instars. 
Although the resolving power of the In- 
transformed data set nearly doubled the 
number of discernible groups, classification 
accuracy remained high at 95.3% (i.e., 
4.7% of cases were assigned to the incorrect 
group). Clearly, the possibility exists that 
additional instars exist that were not dis- 
cerned. Nonetheless, by using a suite of In- 
transformed variables, this method of instar 
determination represents a significant im- 
provement over the commonly used head 
capsule width analysis because both resolv- 
ing power and precision are substantially 
increased. 

Insects with a small number of distinct 
instars may have greater inter-instar vari- 
ance in size than intra-instar variance, 
which allows for the number of instars to 
be effectively discerned using techniques 
such as the multivariate analysis of body 
measurement presented here. For example, 
Sites (1991), using such methods, con- 
firmed that Pelocoris poeyi (Guerin Mene- 
ville) (Heteroptera: Naucoridae) has five 
nymphal instars. In mayflies, overall in- 
crease in size of the nymph occurs at each 
molt as well as a differential growth rate of 
body parts (allometric growth); therefore, 
one function of molting is to change mor- 
phological structure (Ide 1935). Although 
the method discussed here was shown to 
improve our ability to discriminate among 
instars from a field population, laboratory 
rearing of mayflies would provide the most 
accurate data, although environmental influ- 
ences on body size and shape and instar 
number would not be realized. 

Population studies of species of Steno- 
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nema have revealed mainly univoltine win- 
ter cycles, which are characterized by a sin- 
gle generation overwintering in the nymph- 
al stage (Clifford 1982). The pattern ob- 
served during this research for S. 
mediopunctatum approximates this pattern, 
with overwintering nymphs, continued 
growth throughout spring, and emergence 
in early summer. The wide variation in fre- 
quency values and inconsistent patterns 
most likely indicate a large number of 
nymphal instars present in the stream at a 
given time. 

Temperate populations of Tricorythodes 
generally have been characterized as mul- 
tivoltine (Clifford 1982). Although some 
temporal size shifts were observed during 
this study, a pattern was not clear. Because 
multivoltine life history patterns typically 
would be evidenced by a wide range of siz- 
es of individuals throughout the year, T. 
curvatus apparently also is multivoltine in 
this region. 

Life histories of many Caenis popula- 
tions have been reported to be quite flexible 
(Clifford 1982). Nearly half of the Caenis 
species for which life histories have been 
documented are univoltine winter and half 
are bivoltine winter-summer, with an over- 
wintering generation in the nymphal stage 
and one summer generation (Clifford 
1982). For both Caenis species examined 
here, the largest individuals occurred from 
late April through mid-June, with a summer 
emergence period over four months (June— 
September). Therefore, it appears that both 
C. hilaris and C. latipennis have a univol- 
tine life history with an overwintering 
nymphal stage. Black light collections of 
adults of these species also support univol- 
tinism. 

Geographic variation in voltinism due to 
thermal or nutritional regimes commonly 
has been reported (Sweeney 1984); how- 
ever, even within populations, variability in 
life cycle duration can exist (Butler 1984). 
Therefore, because of the potential for en- 
vironmentally induced variation that exists 
for the majority of aquatic insects (Wallace 
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and Merritt 1980), life history attributes, 
such as instar number and voltinism, should 
not necessarily be considered species spe- 
cific traits, but should be examined on a 
regional basis for each population. Larger 
sample sizes may have improved resolution 
in the size-frequency graphs, although a 
high degree of variation is inherent, and 
should be expected, in mayfly populations. 
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